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. Abstract. Three pulsar timing arrays are now producing high quality data sets. As reviewed 

T-H ■ in this paper, these data sets are been processed to 1) develop a pulsar-based time standard, 2) 

' search for errors in the solar system planetary ephemeris and 3) detect gravitational waves. It is 

■ expected that the data sets will significantly improve in the near future by combining existing 
-4^ ■ observations and by using new telescopes. 
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1. Introduction 

Pulsar timing data sets are now of sufficient length and precision to start to realise 
many of the goals of "pulsar timing arrays" (PTAs). The first major PTA was initiated 
" in 2004 using the Parkes radio telescope (Manchester et al., 2012) and is known as the 
7^ Parkes Pulsar Timing Array (PPTA) project. The European PTA (EPTA) makes use 
O ' of the telescopes at Jodrell Bank, Westerbork, Effelsberg, Nangay and Sardinia (e.g., 
^ \ Ferdman et al. 2010). The North American PTA (NANOGrav; Jenet et al. 2009) obtains 
^ ■ observations using the Arecibo and Green Bank radio telescopes. Together these three 
I— —I. PTAs form the International Pulsar Timing Array (IPTA; Hobbs et al. 2010a) which pro- 
vides high quality timing observations of approximately 40 of the most stable millisecond 
" ^ , pulsars known. 

■ The first stage of PTA-related data analyses is to determine the pulse times-of-arrival 
[~--. , (ToAs) for each observation of each pulsar. These ToAs are converted from the observa- 
C ' tory time standard to a realisation of terrestrial time (TT). Barycentric arrival times are 

, calculated using knowledge of the relative position of the Earth with respect to the solar 

■ system barycentre using a planetary ephemeris and by converting from TT to coordinate 
, barycentric time (TCB). These barycentric arrival times are compared with predictions 

of the arrival times using a model for the pulsar rotational and orbital parameters. The 
. . , differences between the actual measurements and the predictions are known as the "pul- 
^ ■ sar timing residuals" . This technique, known as "pulsar timing" , is widely used in pulsar 
astronomy and is described in detail by Edwards, Hobbs & Manchester (2006). 
' Various phenomena such as gravitational waves, unexplained timing irregularities, 

, glitch events, errors in terrestrial time standards or in the solar-system ephemeris will 
induce timing residuals. The main aim of PTAs is to distinguish between these various 
phenomena by searching for correlations between the timing residuals of multiple pulsars. 
For instance, pulsar timing irregularities, glitch events or interstellar medium variations 
will lead to timing residuals that are uncorrelated between different pulsars. In contrast 
an error in the terrestrial time standard will lead to timing residuals that are identi- 
cal for different pulsars (assuming that all pulsars have been observed over the same 
time span). Errors in the solar-system ephemeris will affect different pulsars depending 
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Figure 1. This figure is reproduced from Hobbs et al., (2012). Tlie top panel siiows the sampling 
for each of the pulsars in our sample. The lower panel shows the difference between the pulsar 
timescale and TT(TAl) as points with error bars. The solid line indicates the difference between 
TT(TAI) and TT(BIPMll) after a quadratic polynomial has been fitted and removed. Full 
details are available in Hobbs et al. (2012). 

upon their ecliptic coordinates. The phase and amplitude of timing residuals induced 
by a gravitational wave will depend upon the pulsar-Earth-source angle (e.g., Detweiler 
1979). The expected correlation between different pulsars for timing residuals induced by 
an isotropic, stochastic gravitational wave background has been calculated by Hellings 
& Downs (1983). 



2. Developing a pulsar time standard 

Millisecond pulsar rotation is incredibly stable. This leads to the possibility of develop- 
ing a time scale based on the pulsar rotation analogous to the free atomic scale, Echelle 
Atomique Libre (EAL). The Ensemble Pulsar Scale (EPS) can be used to detect fluc- 
tuations in atomic timescales and therefore can lead to a new realisation of Terrestrial 
Time, TT. 

Earlier attempts to develop a pulsar timescale have been made by Guinot & Petit 
(1991), Petit & TaveUa (1996), Rodin (2008) and Rodin & Chen (2011). Recently we 
have developed a method to produce a new time scale based on observations of 19 pul- 
sars obtained for the PPTA project (Hobbs et al. 2012). The new algorithm has been 
implemented as part of the tempo2 software package (Hobbs, Edwards & Manchester 
2006). This algorithm accounts for various features of the observations such as: 1) irreg- 
ular sampling, 2) different data spans for different pulsars and 3) different fitting param- 
eters for different pulsars. Our result is reproduced in Figure [T] We successfully follow 
features known to affect the frequency of the International Atomic Timescale (TAI) and 
we find marginally significant differences between our pulsar time scale, TT(PPTAll), 
and TT(BIPMll). 

This work is being continued by combining the Parkes observations with data from 
other observatories. The new analysis will confirm or deny the tentative discrepancies 
between TT(PPTAll) and TT(BIPMll) whilst significantly improving the stability and 
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Figure 2. Offset of tlie Eartli-Solar System barycentre vector as a function of time from the JPL 
DE421 planetary epiiemeris. Tlie tliree panels show the offset in the three spatial coordinates. 

precision of the pulsar scale. In the longer term it is expected that a future pulsar time 
scale will be combined with the best atomic timescale to give the world's most stable 
time scale that will be valid effectively forever. 



3. Improving the solar system ephemeris 

The pulsar timing method relies on the determination of pulse times of arrival as mea- 
sured in the solar-system barycentre. The procedure requires knowledge of the position 
of the Earth with respect to the solar-system barycentre. This is obtained using a pub- 
lished solar system ephemeris. Errors in the ephemeris will lead to timing residuals. For 
instance, an error in the mass of Jovian system assumed when forming the ephemeris 
will lead to residuals proportional to both the pulsar-barycentre- Jupiter angle and the 
size of the mass error. As Jupiter orbits the barycentre, the angle will change and hence 
sinusoidal pulsar timing residuals will be induced with a period equal to that of Jupiter's 
orbit. 

Champion et al. (2010) searched for such sinusoidal timing residuals using PTA obser- 
vations of four pulsars obtained using the Arecibo, Parkes and Effelsberg radio telescopes. 
In most cases published masses obtained from space-craft data were more precise than 
the pulsar results. However, for the Jovian system, the Champion et al. (2010) measure- 
ment of 9.547921(2) x 10^"* Mq is more accurate than the mass determined from the 
Pioneer and Voyager space-craft. 

The Champion et al. (2010) technique can only be applied to known solar-system 
objects. However, it is also possible to determine offsets from the predictions of a specific 
ephemeris in the Earth's position with respect to the barycentre. Significant offsets in 
any of the three spatial coordinates can subsequently be analysed in order to identify the 
orbital parameters of any previously unknown object. In Figure [2] we show the results 
of an initial analysis using the PPTA observations. We plot the offset in the Earth- 
barycentre vector as a function of time compared with the value predicted using the JPL 
DE421 Solar System ephemeris. We identify no significant offsets suggesting that the 
ephemeris is adequate for our current purposes over our five year data span. 
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Figure 3. Current upper bounds on the gravitational wave background and past and current 
predictions. The different symbols are shown in the text. 

The sensitivity of a PTA to errors in the solar system ephemeris depends upon the 
timing precision achieved and the data span. It is expected that significant improvements 
will occur when the data span becomes longer than 29 years, the orbital period of Saturn. 
Combining the observations from the existing timing arrays will also significantly improve 
our sensitivity. 

4. Searching for gravitational waves 

Sazhin (1978) and Detweiler (1979) showed that gravitational waves (GWs) passing 
through the solar system will induce timing residuals that are potentially detectable us- 
ing PTAs. Pulsar data sets are sensitive to GWs with periods longer than the typical 
data sampling and shorter than the total time span of the observations. Hence, pulsar 
experiments are sensitive to ultra-low frequency (~ 10~^-10~^ Hz) GWs and are com- 
plementary to ground-based and space-based GW detectors such as LIGO and eLISA. 

Sources of a background of GWs include cosmic strings (see Sanidas et al. 2012, Regim- 
bau et al. 2012 and references therein), the inflationary era (e.g., Zhao 2011) and coa- 
lescing supermassive binary black holes (e.g., Sesana, Vecchio & Colacino 2008). The 
induced timing residuals induced by a GW background is often described as having a red 
power spectrum: 



The almost horizontal lines in Figure [3] show recent upper bounds that have been placed 
on A(ol) by Jenet et al., (2006), Van Haasteren et al. (2011) and Demorest et al. (2012). 
For comparison an earlier upper bound by Kaspi, Taylor & Ryba (1994) is plotted as a 
downward pointing arrow. 

Since the start of the PPTA project our expectations for the detectable GW signal have 
changed. In Figure[3]we show the past and current predictions for the signal strength and 
spectral exponent. Estimates that were available at the time that our first upper bound 
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was published (Jenet et al. 2006) for cosmic strings, the inflationary era and black holes 
are shown as the left-most vertical solid line, the hollow rectangle and the right-most 
vertical solid line respectively. Since 2006, new research has (1) broadened the possible 
range of acw indicated that very few constraints are available on the lower amplitude 
of the cosmic-string background signal (Sanidas, Battye & Stappers 2012; shown as the 
left-most hashed region in Figure [3]), (2) confirmed that the most-likely signal from the 
inflationary era has a very low amplitude (the small hashed region near the bottom-centre 
of the figure should be considered as an upper-bound on the amplitude of this signal) and 
3) led to the possibility that the background caused by coalescing black holes is fiatter 
than originally predicted (right- most hashed region). 

Earlier work has assumed that the detected GW signal will be an isotropic, stochastic 
background. More recent studies have suggested that the exact nature of the detected 
signal could be a background, but may be an individual non-evolving source, a chirping 
system, a memory event, or a burst event. Various algorithms have therefore been devel- 
oped to search for these various signals (e.g., Yardley et al. 2010; Finn & Lommcn 2010; 
Cordes & Jenet 2012). 



5. Improving the data sets 

The sensitivity of a PTA data set to a given science goal depends upon the timing 
precision achieved, the noise present in the data, the number of pulsars observed and 
the data span. For most millisecond pulsars over decadal time scales the dominant noise 
source is caused by turbulence in the interstellar medium leading to variations in the 
pulsar's dispersion measure (e.g.. You et al. 2007). This effect can only be removed by 
observing the pulsar at widely separated observing frequencies. To help address this 
problem many PTAs are now developing wide-band receiver systems that provide a large 
frequency coverage for each observation. 

For many pulsars, the dominant uncorrectable noise is caused by receiver noise and 
pulse jitter. Receiver noise can only be reduced using more sensitive telescopes or longer 
observations. In the future it is expected that telescopes such as the Five-hundred metre 
spherical telescope (FAST) in China, or the Square Kilometre Array (SKA) will provide 
a huge increase in sensitivity. However, individual pulse-shape variations or jitter may 
provide a limit to the precision with which pulse arrival times can be measured and 
therefore require a modified strategy in the use of these new telescopes (Oslowski et al. 
2011, Liu et ah 2012). 

Over long time scales pulsars are known to exhibit irregularities in their spin-down 
rate (e.g. Verbiest et al. 2009 and Hobbs et al. 2010b). It is currently thought that 
this noise is uncorrectable and will limit the stability of pulsars over long time scales. 
However, recent work (Lyne et al. 2011) has shown that it may be possible to identify a 
deterministic component to these irregularities which opens up the possibility of at least 
partially correcting for their effects. In any case, the discovery of new, stable pulsars is 
necessary to improve the sensitivity of the PTA projects. Numerous surveys are ongoing 
(e.g., Cordes et al. 2006, Keith et al. 2010, Boyles et al. 2012) and are leading to the 
discovery of a large number of new, millisecond pulsars. Including these new pulsars in 
existing PTAs and the likelihood of a large number of new, sensitive radio telescopes in 
the relatively near future suggest that the future is bright for PTA research. 
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